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Abstract 



We use lookback time versus red shift data from gala xy clusters (|CaDOzziello et alJ . 120041 ) 
and passively evolving galaxies ([Simon et al. and apply a bayesian prior on the 

total age of the Universe based on WMAP measurements, to constrain dark energy 
cosmological model parameters. Current lookback time data provide interesting and 
moderately restrictive constraints on cosmological parameters. When used jointly with 
current baryon acoustic peak and Type la supernovae apparent magnitude versus redshift 
data, lookback time data tighten the constraints on parameters and favor slightly smaller 
values of the nonrelativistic matter energy density. 

Keywords: cosmological parameters, galaxies: general, galaxies: clusters: general 
PACS: 95.35.+d, 98.80.-k 



1. Introduction 



It is now a well established fact that the expansion of the Universe is accelerating, 
but the underlying mechanism which gives rise to this cosmic acceleration is still a mys- 
tery. Recent cosmol ogical observations i ncluding the Hubble diagram of Type la super - 



novae (SNela, e.g., Ifficken et all l2009t IShafieloo et"all . l2009t ICuimaraes et all l2009l) 



combined with cosm i c m icrowave background (CMB) anisotropy measurements (e.g 
Dunklev et al.l . 120091: iKomatsu et al... 2009), baryon acoustic peak galaxy power spec 



trum data ( e.g., Percival et al.l 2007 : Samushia fc Ratral . l2009al : ICaztafiaga et al 



Allen et al 



2009; 



2008 : 



Wang, 2009), and galaxy cluster gas mass fraction measurements (e.g 
Sam ushia & Ratra, 2008; Et tori et al.. . 2009.) indicate that we live in a spatially-flat uni- 
verse where nonrelativistic matter contributes about 30% of the critical density. Within 
the framework of Einstein's general theory of relativity, the rest of the 70% of the energy 
density of the Universe is termed dark energy, a mysterious component with negative 
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effective pressure that is re sponsible for the observe d accelerated expansion|l| Fo r recen t 
reviews o f dark ener gy see iRatra fc Voeelev ( 2008 ). Caldwell fc Kamionkowskil (j2009ll . 
FriemanI (|2009l) . and|Sanii| (|2009l) . 

There are many dark energy candidates. The simplest is Einstein's cosmological 
constant A. In addition, there are other options like XCDM, a slowly rolling scalar field, 
Chaplygin gas, etc., which can also give rise to an accelerated expansion of the Universe. 
In this paper we constrain the parameters of three different dark energy models. The first 
model is the cosmological constant dominated cold dark matter (ACDM) model ([Peebles, 



19841 ). In this model the energy density of the vacuum (the cosmological constant) does 



not vary with time and it has a negative pressure characterized by = ^PA, where pA 
is the vacuum energy density. 

Secondly, we consider the XCDM parameterization of dark energy. In this case dark 
energy is assumed to be a fluid satisfying the following relation between pressure and the 
energy density, = WxPx, with Wx < 0; this is not a physically complete model. Lastly, 
we study the slowly rolling dark energy scalar field model (0CDM) with an inverse 
power-law potential energy density for the scalar fie l d, V o c 4)~°' where a is a nonnegative 
constant (Peebles fc Ratral . Il988l: iRatra fc Peeblesl . Il988l: IPeebles fc RatraL [2003)0 We 
only consider the spatially-flat (j)CDM and XCDM cases. The (j)CDM model with a — 
and the XCDM model with ij^ — —1 are equivalent to the spatially-fiat ACDM model 
with the same matter density. In all three models the nonrelativistic matter density is 
dominated by cold dark matter. 

In this pap er we use two sets of lo okback time versus redshift meas urements, for 



galax y clusters ([Caoozziello et al.l . 120041 ) and for passively evolving galaxies ([Simon et al 



l2005l) . and apply a bayesian prior on the total age of the Universe based on WMAP esti- 
mates ([Dunklev et all . [2009*) . to constrain parameters of these dark energy models. This 
time-based cosmological test differs from other widely-used distance-based cosmological 
testsd An important feature of this time-based method is that the age of distant ob- 
jects are independent of each other. Therefore, it may avoid biases that are present in 
techniques that use distances of primary or secondary indicators in the cosmic distance 
ladder method. In the literature a variety of time-based methods have been considered, 
based on measurements of the absolute age of objects, differential age of objects, and 
lookback time of objects (f| 



^For discussions of modification of Einsteinian gravit y on cosmological scales t hat attempt to do away 
with t h e need for dark ene r gy, see iRapetti et al.l l l2009h . iBamba fc Gen's | |2009|) . [cTpozziello fc Salzanol 
ll2009l) . rWu fc ChenI ll2009l '). Iztiou et al.l ll2009l '). and references therein. 

^In the (/)CDM model we consider here, (/> only couples gravitation ally to othe r compone nts. For 
model s where 4> also interacts mo r e dire ct ly with other compone n ts, see IChen et al.| ([2 009'). Baldietal 



, iBento fc Gonzalez Felipe ||2009| ). iPettorino et al.l | |2009|) . iGavela et al.l ||2009|) iLa Vacca et al 
and references therein. For othe r dark e nergy models, s ee iLu et al.ll2009l). lArbevI ||2008|) . 



Hrvcvna fc Sz vdlowski' f 20091), iBasilakod ||20090 . iTsuiikawa et al.l ||20090 . I Dutta^ Scherren 



20091) . 



Neupane fc T rowland (2009), and references therein. 



''Distance-based cosmological tests include those mentioned above that use SNela, CMB, baryon 
acoustic peak, and galaxy cluster gas mass fraction data, as we l l as radio-galaxy and quasar angular size 
versu s redshift data (e.g.,|Chcn & Ratra. 2003': Poda riu et all 12003': Daly ct al., 2009; Santos fc Lirnal 
20081) and gamma-ray burst l uminosity distance versu s redshift m easurements (e.g., Tsutsui ct al., 20091: 



Qi et all I2OO9I : IWei fc Zhang . 2008; Liang fc Zhangj . 12008: Wang, 2008; Sa musliia fc Ratral. I2010D . 



•^A variation of this test uses measurements of the Hubble parameter as a function of redsh ift (e.g., 
ISamushia fc Ratral. l2006l : iLin et al [20081 : iDev et al.L l2008l : iFernandez-Martinez fc Verdd. l2008l . and ref- 
erences therein). 
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The absolute age method is based on the simple criterion that the age of the Universe 
at a give n redshift is always greater than or equal to the age of the oldest object at that 
redsh ift (|Alcaniz fc Limal . ll999t iLima fc Alcanij . I2OOOI: Ijain fc Devi . I2OO6I: IWei fc Zhang , 
20071 ). The differential age method is based on the measurement of Az/Af. Az is the 



red shift separation between the two passively ev olving galaxies having the age difference 
At ( Jimenez fc Loeb , 2002 ; Jimenez et al ]. l2003h . This method requires a large sample of 
passively evolving galaxies with high quality spectroscopy and is probably more reliable 

than the absolute age method as a number of systematic effects are elimin ated. 

Lo okback time as a tool to constrain dark energy models was first used bv lCapozziello et al 



(|2004[ ) who compiled a list of galaxy cluster ages and redshifts and used this data to con- 
strain the XCDM dark energy parameterization. This data has been used to constrain 
brane cosmology and holographic dark energy models ()Pires et al.L I2OO6I : lYi fc Zhanj . 



20071 ). The lookback time test has also been applied using passively evolv ing galaxies 
data, to constrain parameters of XCDM and ACDM (jDantas et al.l . l2007l [2009). No 
doubt these time-based methods are subject to some different systematic errors but they 
offer an independent means to cross-check cosmological constraints obtained using other 
techniques. 

In this paper we t ake advantage of th e fact that the Capozziello et al. ( 20041) galaxy 
cluster data and the Simon et al. I (l2005l) passive galaxy data are independent, so it is 
straightforward to use them simultaneously in a lookback time versus redshift test anal- 
ysis of dark energy models. Our joint analyses of these data sets allow us to derive the 
tightest lookback time constraints on dark energy parameters to date. The resulting 
constraints are moderately restrictive, and these data favor lower matter density values 
than do some other current data, but are consistent with a spatially-flat ACDM model 
in which nonrelativistic matter contributes 30% of the energy budget at a little less than 
two standard deviations. To derive tighter constraints, we perform a joint analysis of the 
lookback time data with current baryon acoustic peak and SNela measurements. 

In Sec. 2 we describe the lookback time as a function of redshift test. The data and 
method we use are outlined in Sec. 3. Our results are presented and discussed in Sec. 4. 



2. Lookback time versus redshift test 



The lookback time is the difference between the present age of the Universe (tg) and 
its age at redshift z, t(z), 



tL{z,p) 



Hp) - - ^ 

no 



dz' 



dz' 



{1 + z')n{z',p) 



dz' 



{l + z')H{z',p) 



Ho Jo {l + z')H{z',py 



(1) 



Here p are the parameters of the cosmological model under consideration, 'H{z^p) = 
H{z,p)/ Hq, H{z,p) is the Hubble parameter at redshift z, and the Hubble constant 

Ho = lOO/i k m s~^Mpc~^ 

Following ICapozziello et ah ( 2004 ). the observed lookback time t°^^{zi)^ to an object 
i at redshift z,- is defined as 



V^4; ''inci ''0 ) — ''0 ^i\Zi) ''1] 
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(2) 



Here 

• tQ^^ is the measured current age of the Universe. 

• ti{zi) is the age of the object (passively evolving galaxy, cluster, etc.), defined as the 
difference between the current age of the Universe at redshift Zi and the age of the 
Universe when the object was born at redshift Zf, 

1 dz' 

U{z.) = t{z.) - t(z,) = t,[z,) t,iz.) = ^ J^^ (l + ,OH(z',rt ' 

where we have used Eq. (1). 

• tine — tQ^^ — tL{zf) is the incubation time of the object. This delay factor encodes our 
ignorance of the formation redshift Zf. 

To compute model predictions for the lookback time tL{z,p), Eq. (1), we need an 
expression for H{z,p). In the ACDM model the Hubble parameter is 

Hiz,p) = Ho [a„(l + zf + (1 - - 1]a)(1 + z)2 + Qa]^/^, (4) 

where p are r2,n and JIa, the nonrelativistic matter and dark energy density parameters 
at z = 0. For the XCDM parameterization in a spatially-flat cosmological model we have 

H{z,p) - Ho [f^n,(l + zf + (1 - a„)(l + z)3(i+-J]i/2, (5) 

where p are fim and Wx- In the spatially- flat (/>CDM model 

Hiz,p) = Ho [r!„,(l + zf + n^z)]^/^, (6) 

where the scalar field energy density parameter fltf,{z) can be evaluated numerically by 
solving the coupled set of equations of motion, 

+ 3^0-^m>-("+i)=O, (7) 
a 2. ^ 

= ^[^r.{l+zf + n^{z% (8) 

n^{z) = + Kmlr'']n2. (9) 

Here a{t) is the scale factor, an overdot denotes a time derivative, nip is Planck's mass, 
and K and a are non-negative constants that characterize the inverse power law potential 
energy density of the scalar field, V{(j>) = k(I)~°'. In this case the parameters p are fim 
and a. 



3. Data and computation 



In order to constrain cosmologic al parameter s of AC DM, XCDM, and cjjCDM, we 
use two age data sets. One is the ISimon et al. ages of 32 passively evolving 

galaxies (Table 1, R. Jimenez, private communication 2007) in the redshift interval 
0.117 < z < 1.845. For this sample we assume a 12% one standard deviation un- 
certainty on the age measurements (R. Jimenez, private communication 2007). The 
other is the ICapozziello et al. ()2004l. Table 1) ages of 6 galaxy clusters in the redshift 
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range 0.10 < z < 1.27. This sample has a 1 Gyr one standard deviation uncertainty on 
the age measurements. In all, we have 38 measurements of ^^^^{zi) with uncorrelated 
uncertainties tXi. 

For each model and parameter value set (p) we compute the function 



38 

E 

=1 



(10) 



where a°^^ is the uncertainty in the estimate of to and tL[zi,p) and tQ{p) are the predicted 
values in the model under consideration. From construct a likelihood function 

£'(p,ifo,tinc) cxexp(-xV2). 

The likelihood function i/o, tine, ^o^*^) depends on the total age of the Universe 
, incubation time tine and the Hubble parameter Hq. We do not kno w t\-nr. and so treat 



it as a nuisance parameter a nd analytically marginalize C over it as in lCapozziello et al 



(j2004r ) ; iDantas et al.l pOOTI ) . We treat Hq as a n uisance parameter and marginalize over 
it with a Gaussian prior with h = 0.742 ± 0.036 ( Riess et al. . 2009f) . a little higher than 
but s till consistent with, the earlier summary value oi h = 0.68 ± 0.04 (jGhen et al 



20031). We also also apply a bayesian prior as a Gaussian function with central values 
and variances based on the WMAP estimate of the total age of the Universe, which 
is tf"" = (13.75 ± 0.13) Gyr for the A CDM model and t°'°'' = (13.75l^;2?) Gyr for 
the XCDM model (jPunklev et al.l . |2009|) FI For the ^CDM model we assume the same 
central value as the other two models and conservatively inflate the error bar to t^^^ = 
(13.75 ± 0.5) Gyr. The resulting lookback time likelihood function depends only on 
the two cosmological parameters p, Cl{p)- The best fit parameters are the pair p* that 
maximize the likelihood function and the 1,2, and 3f7 confidence level contours are defined 
as the sets of cosmological parameters Pa at which the likelihood C{pa) is exp(— 2.30/2), 
exp(— 6.18/2), and exp(— 11.83/2) tim es smaller than th e maximum likelihood C{p*). 

To check our method we used the ICapozziello et al.l ()2004h galaxy cluster data and 
the earlier result they used and computed the constraints on th e XGDM parame- 
teriza tion. Our contours are consistent with those shown in Fig. 2 of ICapozziello et~al 



(12004). We also used thelSimon et al.l ()2005[ ) passively evolving galaxy ages and the 
value iDantas et al.l (|2007l ) used to constrain the ACDM model. We find that if we pick 
h = 0.72 we are ab le to accurately reproduce the central and right panels of Fig. 2 of 
Dantas et all (|2007l ). 

The lookback time versus redshift data constraints on ACDM, XCDM, and (/)CDM 
are shown in Figs. 1-3. 



4. Results and discussions 

Figure 1 shows the constraints on the ACDM model from the lookback time and age 
of the Universe measurements. The data favor low vales of both fim and JIa with the 
best-fit values being ilm = 0.01 and JIa = 0.19. These data prefer spatially-open models, 
however a spatially-flat ACDM model with fim — 0.3 is less than 3cr from the best fit 
model. The data constrains fim to be less than 0.45 on ia confidence level. 



^The numbers are taken from [http:/ /lambda. gsfc.nasa.gov/[ 
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Figure 2 presents the constraints on the XCDM parametrization of the equation of 
state. The nonrelativistic matter density parameter is constrained to be less than 0.5 at 
3(7 confidence. Low values of rim are favored with the best-fit values being = 0.03 
and = —0.41 and a spatially-flat model with Jim — 0.3 is about 2(t from the best fit 
model. 

Figure 3 shows the constraints on the 0CDM model of dark energy. In this model 
the nonrelativistic matter density parameter is less than 0.5 at 3a confidence. The a 
parameter on the other hand is not well constrained. The best-fit parameter value is 
a = 10, but the likelihood is very flat in the direction of a and the difference between 
the best-fit value and a = (which is the spatially-flat ACDM case) is slightly less than 
2a. 

Current lookback time data by themselves are unable to tightly constrain cosmological 
parameters. C o nstra ints from galaxy cluster gas mass fraction versus redshift data (e.g., 
Chen fc Ratra . 2004i) . SNela apparent magnitude versus r edshift measuremen t s (e.g 



Wilson et al.l 



2006 ). and baryon acoustic peak data (e.g., Samushia fc Ratra . l2009b ) 
are more restrictive than the lookback time constraints. However, the constraints from 
lookback time dat a are somewhat ti ghter than the constraints from strong gravitational 

2004. parameter as a func- 

20071) . radio galaxy angular size versus redshift 



lensing data fe.g.. iChae et al, 
tion of redshift (e.g., 



. , .Samushia et al.l „ ^ 

data (e.g., Dalv et al.l, l2009l), and g amma-ray burst luminosity distance versus redshift 
data (e.g.. ISamushia fc Ratra . 2010). 

To get tighter constraints on cosmological parameters we combine the lookback time 
data and the measur ement of the age of the Universe with baryon acoustic peak data 
(IPercival et al. . 2007 ) and SN ela "Union" apparent magnitude versus redshift measure- 



ments iKowal ski et al. . 2008f) . Since these data sets are independent we compute a joint 



likelihood function that is a product of individual likelihood functions 



SNc, 



(11) 



and define the best-fit parameters and confidence level contours as discussed above. 

The constraints on the three dark energy models from a joint analysis of these data are 
shown in Figs. 4-6. Currently available lookback time data do not significantly change 
the results derived using BAO peak measurements and SNela apparent magnitude data. 
In all three dark energy models when lookback time data are added to the mix the 
confidence level regions favor slightly smaller values of nonrelativistic matter density 
parameter 

Overall, current data is a good fit to all three dark energy models. For (/)CDM 
and XCDM they slightly favor time-dependent dark energy, but the time-independent 
cosmological constant is also a good fit. 

We anticipate that a new, improved data set of lookback times will soon be available 
(R. Jimenez, private communication, 2009.) With more and better data we expect 
significantly tighter constraints on dark energy parameters. The lookback time versus 
redshift test, either by itself or at least in combination with other cosmological probes, 
could prove very useful in detecting or constraining dark energy time evolution. 
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Figure 2: 1, 2, and 3(7 confidence level contours for the XCDM parameterization of dark energy in 
a spatially-flat cosmological model, from the lookback time data and measurement of the age of the 
Universe. The dashed u)x = — 1 line corresponds to spatially-flat ACDM models. The cross indicates the 
best-fit parameters fJm = 0.03 and oJx = —0.41 with = 28 for 37 degrees of freedom. 
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Figure 3: 1, 2, and 3(7 confidence level contours for the spatially-flat </)CDM model from the lookback 
time data and measurement of the age of the Universe. The a = horizontal axis corresponds to 
spatially-flat ACDM models. The cross indicates the best-fit parameters Hm = 0.04 and a = 10 with 
= 22 for 37 degrees of freedom. 
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Figure 4: 1, 2. and 3cr confidence level contours for the ACDM model. Numerical noise is responsible 
for the jaggeduess of parts of the contours. The dashed line demarcates spatially-flat models. Dotted 
lines (circle denotes the best-fit point at Qm = 0.30 and Q/^ = 0.78 with = 359 for 346 degrees of 
freedom) are derived using the lookback time data, measurement of the age of the Universe, SNela Union 
data, and BAO peak measurements, while solid lines (cross denotes the best-fit point at Om = 0.32 and 
Ha = 0.78 with = 318 for 307 degrees of freedom) are derived using SNela and BAO data. The 
dashed line corresponds to spatially-flat models. 
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Figure 5: 1, 2, and 3<t confidence level contours for the XCDM parameterization of dark energy in a 
spatially-fiat cosmological model. The dashed line demarcates spatially-flat ACDM models. Dotted 
lines (circle denotes the best-fit point at Qm = 0.19 and lo^ = —0.80 with = 352 for 346 degrees of 
freedom) are derived using the lookback time data, measurement of the age of the Universe, SNcIa Union 
data, and BAO peak measurements, while solid lines (cross denotes the best-fit point at Hm = 0.19 and 
= —0.81 with = 321 for 307 degrees of freedom) are derived using only SNela and BAO data. 
The dashed Wx line corresponds to spatially-flat ACDM models. 
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Figure 6: 1, 2, and 3a confidence level contours for the spatially-flat <j>CT>M model. The a = horizontal 
axis corresponds to spatially-flat ACDM models. Dotted lines (circle denotes the best-flt point at f2m = 
0.215 and a = 0.0 with = 359 for 346 degrees of freedom) axe derived using the lookba<;k time data, 
measurement of the age of the Universe, SNela Union data, and BAO peaJc measurements, while solid 
lines (cross denotes the best-fit point at Clm = 0.22 and a = 0.0 with = 329 for 307 degrees of 
freedom) axe derived using only SNela and BAO data. 



14 



Simon et al. (2005) Ealaxv aees 


Zi 


U{z^) (Gyr) 


0.1171 


10.2 


0.1174 


10.0 


0.2220 


9.0 


0.2311 


9.0 


0.3559 


7.6 


0.4520 


6.8 


0.5750 


7.0 


0.6440 


6.0 


0.6760 


6.0 


0.8330 


6.0 


0.8360 


5.8 


0.9220 


5.5 


1.179 


4.6 


1.222 


3.5 


1.224 


4.3 


1.225 


3.5 


1.226 


3.5 


1.340 


3.4 


1.380 


3.5 


1.383 


3.5 


1.396 


3.6 


1.430 


3.2 


1.450 


3.2 


1.488 


3.0 


1.490 


3.6 


1.493 


3.2 


1.510 


2.8 


1.550 


3.0 


1.576 


2.5 


1.642 


3.0 


1.725 


2.6 


1.845 


2.5 



